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^SJ , Abstract 

We investigate the resonant states for the parabolic potential barrier known also as 
inverted or reversed oscillator. They correspond to the poles of meromorphic continuation 
(•f^ , of the resolvent operator to the complex energy plane. As a byproduct we establish an 

' interesting relation between parabolic cylinder functions (representing energy eigenfunc- 

tions of our system) and a class of Gel'fand distributions used in our recent paper. 

> 
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1 Introduction 

In a recent paper [1] we have investigated a quantization of the simple damped system^ 



' u = -'ju . (1.1) 



To quantize this system we double the number of degrees of freedom, i.e. together with (1.1) 
we consider v = +'yv. The enlarged system is a Hamiltonian one and its quantization leads 
to the following quantum Hamiltonian: 

H = -^{mJ + ^^l) . (1.2) 

We showed that the above system displays two families of generalized eigenvectors corre- 
sponding to purely imaginary eigenvalues Hfl^ = ±Enfn- These eigenvectors are interpreted 
as resonant states — they correspond to the poles of energy eigenfunctions when continued to 
the complex energy plane. It turns out that resonant states are responsible for the irreversible 
behavior. We showed that there are two dense subspaces ^± G L^(R) such that restriction 
of the unitary group U{t) = e~*^* to ^± does no longer define a group but gives rise to two 
semigroups: U-{t) = [/(t)|$_ defined for t > and U+{t) = [/(t)|$_^ defined for t < 0. In 

^We slightly change the notation: the coordinates {x,p) used in [1] are replaced by {u,v) in the present 
paper. 
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the framework of Gel'fand triplets (see e.g. [2]) it means that the quantum version of the 
damped system (1.1) corresponds to the Gel'fand triplet: 

C l2(M) C , (1.3) 

together with the Hamiltonian H\^_ . This system serves as a simple example of Arno Bohm 

theory of resonances [3] (see also [4, 5]) and illustrates mathematical results of [6]. 

In the present paper we continue to study this system but in a different representation. 
Let us observe that performing the linear canonical transformation (u,v) — > ix,p): 

-fx-p -fx + p 

u = — , V = — , (1.4) 

one obtains for the Hamiltonian 

H = lif-j'x'). (1.5) 

It represents the parabolic potential barrier V(x) = —^'^x^/2 and it was studied by several 
authors in various contexts [7, 8, 9, 10, 11, 12, 13]. It is well known that this system gives 
rise to the generalized complex eigenvalues — the physical reason for that is the potential 
unbounded from below. We find the corresponding energy eigenstates for (1.5). They are 
given in terms of parabolic cylinder functions Dv{x). Using the Gel'fand-Maurin spectral 
decomposition we find the resolvent operator R{z,H) = {H — z)"^ and relate its poles to 
the resonant states. As a byproduct we established a deep relation between the Gel'fand 
distributions [14, 15] (used in [1]) and parabolic cylinder functions Di,{x). The details are 
included in the Appendix. 



2 Inverted oscillator and complex eigenvalues 

Let us note that H defined in (1.5) corresponds to the Hamiltonian of the harmonic oscillator 
with purely imaginary frequency ui = (in the literature it is also called an inverted or 
reversed oscillator). The connection with a harmonic oscillator may be established by the 
following scaling operator [16]: 

Va := exp Q(xp + px)^ , (2.1) 

with A G R. Using commutation relation [x,p\ = i, this operator may be rewritten as follows 

= e^^P = e-'t g"*^^^- , (2.2) 

and therefore it defines a complex dilation, i.e. the action of oii ^ function (p = ip{x) is 
given by 

Vx (p{x) = e"'^ (fie"'^ x) . (2.3) 

In particular one easily finds: 

VxxV^^ = e-'^x , VxpV^^ = e'^p , (2.4) 



2 



Fa if - I'x') V,-' = e'^' (f - e-^'^Wx^) . (2.5) 
Therefore, for e^*'^ = —1, i.e. A = ±7r/4, one has 

F±,/4^Ft-% = ±z#ho, (2.6) 

where 

^ho = ^(p' + 7'^') , (2.7) 
stands for the oscillator Hamiltonian. In particular if E)^ = 7(71+!) is an oscillator spectrum 

H^^i^l- = , (2.8) 

then 

= ±Enft , (2.9) 

with 

En = iE^'' = i7(^n + ^^ , (2.10) 

and 

ftix) = %^/,i;l^°{x) = e±^f ^P^^ie^'^^x) . (2.11) 
Now, recalling that (see e.g. [17]) 

iP^%x) = e-"2-'Hni^x) , (2.12) 

where stands for the n-th Hermite polynomial and the normalization constant 

~" = (^)'' ''-''^ 

one obtains the following formulae for the generalized eigenvectors of H: 

ftix) = iV± e^^i-'i/„(v^x) , (2.14) 

with 

Clearly, are not elements from (M) but they do belong to the dual of the Schwartz space 
»S(Rj;)', i.e. they are tempered distributions. 

Proposition 1 Two families of generalized eigenvectors satisfy the following properties: 
1. they are conjugated to each other: 



f+(x) = f-(x) , (2.16) 
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2. they are orthonormal 

{fn\fl)=Snm , (2.17) 

3. they are complete 



^f^{x)f^ix') = S{x-x') . (2.18) 

n=0 

The proof follows immediately from orthonormality and completness of oscillator eigenfunc- 
tions V'n"- Formula (2.16) implies that f+ and are related by the time reversal operator 
T: Tijj := ijj. Recall [1] that in n-representation T is unitary (it is defined by the Fourier 
transformation), whereas in x-representation it is antiunitary. 

3 Change of representation 

It should be clear that there exists relation between generalized eigenvectors fn{x) and /^(u) 
found in [1]: 

f+{u) ~ n" , f-{u) ~ ^(")(.i) . (3.1) 

They define the same eigenvectors | ± n ) but in different representations: 

fti^) = {^l^n) , f^{u) = {u\±n) . 

To find this relation let us observe that the canonical transformation (1.4) is generated by 
the following generating function 

S{x, u) = ^x^ - ^xu + , (3.2) 

that is, 

V=— v = -— (3 3) 

9x ' du ' 



Let us define a unitary operator 

by 



U : lH^u) — > 



^Xl 1 



/ {Uf ){x) = C r /(^z)e'^(-'«) du , (3.4) 

J —oo 

where the constant 'C is determined by 

/oo 
-oo 

It implies C = e*" Co, where a is an arbitrary phase and 

C„=(J,)', (3.6) 

In the next section it would be clear that a natural choice for the phase is a = — 7r/8. Clearly, 
U may be extended to act on 5(M„)'. It is easy to show that 

UiSiRuY) C S{R^y . (3.7) 
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Proposition 2 The generalized eigenvectors € 5(M^)' and € 5(Mu)' are related by: 

ft=i^fn- (3.8) 

Proof. Let us show that f+ = W/^, that is 

fn(^) ~ J ^x"e^^(^'")du . (3.9) 
Using the definition of S{x, u) one has 

Now, due to the well known formula for the Hermite polynomials 

e'^'^e--' = (-l)"i7„(z), (3.11) 

one obtains 

' e-^^^lf„(^x) ~ f+(x) . (3.12) 

To prove that f~ = i^/^, let us note that^ 

f;(x) = ~ j u"e-^^(^'")dn . (3.13) 

Now, taking into account that /+ and /~ are related by the Fourier transformation 

= V2^{-i)'^F-^[5^^\k)]{u) , (3.14) 

one obtains 

J ^n^-iS{x,u)j^^ = V2^{-iy' j 5^'''\u)F-^[e-^^]{u)du . (3.15) 

Finally, 

F-^ [e-'^] {u) = / e-'^"" e-'^^'^^^^k = ^e^^(^'") , (3.16) 
v27r J 

and hence 

f;(x) ~ y"(5(")(u)e^^(^'")d« , (3.17) 

which ends the proof. □ 
^It turns out that a function 

S{x,v) = -S{x,v) = -^x'^ + \f^xv - ^v^ , 

serves as a generating function for the canonical transformation (1.4): 

_dS_ _ _dS_ 

dx ' ^ dv 
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4 Energy eigenstates 

The spectrum of the self-adjoint operator (1.2) reads (t{H) = (—00, 00) and the corresponding 
energy eigenstates (in n-representation) are given by (cf. section 6 in [1]): 

with G M. For the basic properties of the tempered distributions G »S(Mu)' we refer 
the reader to [14, 15] (see also the Appendix in [1]). Now, using {x,p) coordinates the 
corresponding eigenvalue problem ^{jp — 7^x^)x^ = Ex^ reads 

dlx^'ix) + (7'^' + 2E)x''{x) = . (4.2) 

Introducing a new variable 

z = ^/2irfx , (4.3) 
the above equation may be rewritten as follows 

5M+(- + ^-^)x^ = 0, (4.4) 

with 

which is the defining equation for the parabolic cylinder functions [18, 19, 20]. Its solution 
X^{z) is a linear combination of D^{z), Dy{—z), D^^^i{iz) and D^^^i(—iz).^ On the other 
hand the energy eigenstates in x-representation x^(^) may be obtained by applying the 
operator U defined in (3.4) to the corresponding eigenstates in u-representation il;^{u): 

/oo 
V'±(^i)e^^(=^'")du . (4.6) 
-00 

Hence 

poo 



e 4 



POO 

/ ee^^-^"/^d^, (4.7) 
Jo 



with y = ^/^2z7x, and using an integral representation for Dp{y) (formula 9.241(2) in [18]):^ 



D.iv) = ^ I r^-' e-y^-i'/' , (4.8) 



^These four functions are linearly dependent. For the linear relation see e.g. formula 9.248 in [18]. 
■'The validity of this formula is restricted in [18] for Rep < 0. However, as we shall show (see the proof of 
Proposition 4), it is valid for all p € C 
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one finds 

Xlix) = Vr^^' T{v + , (4.9) 

with V given in (4.5). Similarly, using an obvious relation {—u)\ = u^, one obtains: 

X^{x) = Vi r{u + l)Z)_,_i(y^x) , (4.10) 

that is, X-{x) = X+{~x)- Actually, instead of x± one may use energy eigenstates with the 
definite parity: 

Xfven = ^(Xf + X-) , (4.11) 

Xfdd = ^(X+-X-) , (4.12) 

that is, 

Xeven Xeven ' Xodd Xodd ' 

(4.13) 

where P stands for the parity operator. 
Proposition 3 Energy eigenstates x± satisfy: 

/oo 

X^{x)xT{x)dx = 6{E-E') , (4.14) 
-oo 

and 

/"CXD 

X±{x)x±ix') dE = d{x - x') . (4.15) 



/ 



The proof follows immediately from the analogous properties satisfied by energy eigenstates 
ijj^ in n-representation [1]. 

In [1] we have used also another generalized basis F[V'i^](u). Now, we find its U image 
in 5(Ma;)'. Recalling the Fourier transformation of (see [14] and Appendix in [1]): 

F[xi]{u) = ^ e^^^i r(A + l){u + iO)-^-^ . (4.16) 
V 27r 

one has 

Therefore, the corresponding x-representation 

nfix) = {U F[^l;-^]){x) , (4.18) 

is given by 

= ^ ^ i2Vir' r(-.) fj^ + ioy e-^^'-^^y^ , (4.19) 
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with y = ■s/2i'yx. Now, using the following integral representation (formula 9.241(1) in [18]) 



D,{y) = —2''+H-ire'^ / + iOy e-^^ 

V7I" 7-0O 



(4.20) 



one obtains 



V^{x) = ' ^ Ti-u)D,{-^2ijx) . (4.21) 



Similarly one shows that 
is given by 

Let us note, that 
and 



r?^(x) = (WF[V'Z^])(x), (4.22) 



r,^{x) = Vi T{-u)D,{^x) . (4.23) 



iy + l = -u , (4.24) 



Vi''^''=Vi . (4.25) 

Clearly, the transition u+l — > —v is equivalent to E — > —E and it corresponds to the fact 
that Hr]^ = —Erj^ while Hx^ = +Ex^- The symmetry between x± ^^nd rj^ fully justifies 
the specific choice of the phase factor in the constant C. One has 



V^{x)=xf{x), (4.26) 

that is they are related by the time reversal operator T: r]±. =T x±- Thus energy eigenstates 
ri^ correspond to the time reversed system. This way all four solutions of (4.4) were used to 
construct four families of energy eigenstates: x+^ X-i V+ ^^'^ V-- 



5 Analytic continuation, resolvent and resonances 

Now, let us continue the energy eigenfunctions x± ^■iid 77^ into the energy complex plane 
E e C and let us study its analyticity as functions of E. 

Proposition 4 The parabolic cylinder function Dx{z) is an analytic function o/A G C. 

For the proof see the Appendix. Due to the above proposition the analytic properties of 
the energy eigenfunctions are entirely governed by the analytic properties of the T function 
which is present in the definition of x± Vi:- Since r(A) has simple poles at A = —n, 
with n = 0,1, 2,..., functions x± have poles at E = —En, whereas functions ri±. have poles 
at E = En, where En is defined in (2.10). Using a well known formula for a residue of the F 
function 

Res(F(A);A = -n) = -^^^ , (5.1) 

n! 

one has 

Res [x^ix); -En) = ^ Vi """^ DniT^/^x) , (5.2) 
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and 

Res (r?|(x); = ^ Vi v^x) . (5.3) 
Hence, using the relation [18, 19, 20]:^ 

L>„(z)=2-te-^iJ„(^-^^ , n = 0,1,2,... , (5.4) 

together with 

Hn{-z) = i-lTHniz) , (5.5) 

one obtains 

Res{xf{xy,-En) ~ ftix) , (5.6) 

and 

Res {ri^{xy,+En) ~ f-{x) . (5.7) 

Now, it is natural to introduce two Hardy classes of functions [21]. Recall, that a smooth 
function / = f{E) is in the Hardy class from above Ti^ (from below H^) if f{E) is a boundary 

value of an analytic function in the upper, i.e. Im£^ > (lower, i.e. Im£^ < 0) half complex 
£^-plane vanishing faster than any power of E at the upper (lower) semi-circle — > oo. 
Define 

:= {4> G 5(M,) I fiE) := {xl\4>) eUl} , (5.8) 

and 

:= [<P e S{R,) I f{E) := (r/||<^) eUl} . (5.9) 

It is evident from (4.26) that = that is 

= T($_) . (5.10) 

Due to the Gel'fand-Maurin spectral theorem [22, 23] any function 0~ G may be decom- 
posed with respect to x± family 

/oo 

dEx^{x){x^\r) , (5.11) 

in 

and any function 0+ G may be decomposed with respect to r/^ family 

/oo 

d£;r?f(x)(r?f|0+) . (5.12) 

in '^^ 

Applying the Residue Theorem one easily proves the following 
^In [18] the corresponding equation 9.253 has a wrong sign. 
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Theorem 1 For any function (f)^ G ^± one has 

oo 

r{x) = ^fn{x){fi\^-) , (5.13) 



n=0 

and 



-(x) = X;e(^)(fnl</>+) ■ (5.14) 



n=0 



The proof goes along the same hnes as the corresponding proof of Theorem 2 in [1]. The 
above theorem imphes the following spectral resolutions of the Hamiltonian: 

/oo 
dEE\xl){xl\ = -J2^n\fn){fi\, (5.15) 
=n n=0 

on and 

/oo °° 
dEE\r,^){r)^\ = J2En\f-}^){a , (5.16) 

I n=0 

on $_|_. The same techniques may be applied for the resolvent operator 

R{z,H) = J—. (5.17) 
H — z 



One obtains 



^Ixf )(x^l = E rF^lfn)(f;ll , (5-18) 



± n=0 

on and 

/oo °° -1 

^l^f = E , (5.19) 

on <I>_|_. Hence, ii(z,//)|$_ has poles at z = —En, and R{z, H)\^^ has poles at z = 

As usual eigenvectors f~ and corresponding to poles of the resolvent are interpreted as 

resonant states. Note, that 

Ri^^H)dz = \ft){f-\:=Pn, (5.20) 
where 7„ is a closed curve that encircles the singularity z = En. Clearly, 

Pn ' Prn — ^nmPn : (5.21) 

and the spectral decomposition of H may be written as follows: 

oo OO 

H = J2 EnPn = - E ^^^n ■ (5-22) 
n=0 n=0 
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Finally, let us note, that restriction of the unitary group U{t) = e defined on the 
Hilbert space L^(M) to no longer defines a group. It gives rise to two semigroups: 



and 



U-{t) : 



U+{t) : $H 



+ ) 



for t>0 , 



for t<0 . 



Using (5.15), (5.16) and the formula for = ij{n + ^) one finds: 

oo 
n=0 

for t>0, and 

oo 



(5.23) 
(5.24) 

(5.25) 



(5.26) 



n=0 



for t < 0. We stress that (/>^ {(t)f ) does belong to L^(M) also for i < (t > 0). However, 
G ((?!)^ G $+) only for t > (t < 0). This way the irreversibility enters the dynamics 
of the reversed oscillator by restricting it to the dense subspace $± of (R) . 

6 Scattering vs. resonant states 

To compare the physical properties of energy eigenstates x± and rj^ and resonant states 
let us investigate its asymptotic behavior at x — > ±oo. Following [20] (see also [10, 11]) one 
finds^ 

X-{x +oo) ~ a/ - exp 



zUx2 + -log(V2^x) + -- + - 



(6.1) 



and 

X^(x -oo) 



X 

— TT — 

— e T 



1 + e" 



-27r- 



exp 



-M^x2 + -log(V^x)--- + - + 



\ - 7 log( V2^x) - 5f + i)] } > (6-2) 



where = argr(— z-^ + 5) ~ ■'^(^ + -'-)• Hence energy eigenstates X- represent scattering 



states (see [10] for more details). The same is true for x+ and 77^. In particular one finds for 



^Putting o = — -E/7 in equation 19.17.9 in [20] and using relation 19.3.1 

V(a,x) = D_^_i{x) , 



one finds: 
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the reflection and transmission amplitudes R and T for x± scattering states [8, 10]: 



£^ 1 m^fl E 



Clearly, computing R and T for time-reversed r]^ scattering states one finds: 



R{rjf) = i?(xf ) , r(r?|) = r(x|) . (6.5) 

Note, that R{x±) ^^'^ ^(x±) have poles at = —En, whereas R{r]^) and T{rj^) have poles 
aX E = +En- Obviously, the corresponding reflection and transition coefficients \R\^ and |rp 
arc time-reversal invariant. 

On the other hand the eigenstates behave as follows: 

f+{x ^ ±oo) ~ (±V^x)"e-^^^' , (6.6) 

and 

f-{x ±oo) ~ (±v/^2;)" e^^^' . (6.7) 

Note, that f~ are purely outgoing states, whereas f+ are purely ingoing states. Moreover, 
resonant states have Breit-Wigner energy distribution. Indeed, 

/CXD 
D,i^ix)f+{x)dx . (6.8) 
-oo 

Now, Di, is an entire function of u and r(— z/) has poles at u = k e N. In the domain where 
n-|-l>Rez^>l one has 



r(— z/) = analytical part -|- } rx 

^-^ k\{k — u)'^ 



(6.9) 



k=0 

Hence, 



(X^|f+) ~ analytical function of ^ + ^ ^ ^ , n ( ffc Ifn ) 

k=o k\ (fc + if + 

~ analytical function of £^ -|- — , (6.10) 

E — En 

which is consistent with the Breit-Wigner formula. 

A Appendix 

The integral formula 9.241(2) in [18] 

/•OO 

-A-1 ^-yi-e/-i 



"rVA) L 



Dx{y) = TTT^ / , (A.l) 
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contains two objects: r(— A) and a distribution ^ ^ ^ which are singular for A = 0, 1, 2, ... . 
However, it is easy to see [14] that 



= (5(-)(0, (A.2) 

p=n 



r(-A) 

which shows that (A.l) defines an entire function of A G C. The same is true for 



oo 



due to 



r(-A) 



= (-ir<5(")(0 . (A.4) 



\=n 

The second integral representation given by 9.241(1) in [18] 

Dx{y) = 2^+2 {-ife^ / + i^f 6-2^ +2^2^^ d^, , (A.5) 

J-00 

where + iO)^ = C+ + e*'^'^^^, seems to have poles at A = —1, —2, .... However, the limit 
limA^_„(^ + iO)-^ is weU defined [14] 

+ ior^ = - -j^zTjT ^^""'^ (0 ■ (A.6) 

Thus, formula (A.5) also defines an entire function of A. 
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